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Abstract 


The symmetry energy at the time of the production of intermediate mass fragments (IMFs) is studied using 


experimentally observed IMF multiplicities combined with quantum statistical model calculations (QSM of Hahn and 


Stócker). The ratios of difference in chemical potentials between neutrons and protons relative to the temperature, 


(u, — Lp) /T , and the double ratio temperature, T, were extracted experimentally in the reactions of 647070, 
94Nj4356N,. 1224s, 1 Au, 23?Th at 404 MeV. The extracted ( H,— Lp) /T scales linearly with dyy, where dyy is the 


asymmetry parameter, (N—Z)/A, of the emitting source and (w, — 4,)/ T =(11.1£1.4)6,, —0.21 was derived. The 


experimentally extracted (4, — 44,)/ T and the double ratio temperatures are compared with those from the QSM 


calculations. The temperatures, T, and densities, p, extracted from the (1, —4,)/T values agreed with those from 


the double ratio thermometer which used the yield ratios of d, t, h and a particles. However the two analyses of the 


differential chemical potential analysis and the initial temperature analysis end up almost identical relation between T 


and p. T=5.25+0.75 MeV is evaluated from the (4, —4,)/T analysis, but no density determination was possible. 


From the extracted T value, the symmetry energy coefficient Esym=14.6+3.5 MeV is determined for the emitting 


source of 7=5.25+0.75 MeV. 
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1 Introduction 


The 


fundamental ingredient in the investigation of nuclear 
[12] 


symmetry energy of nuclear matter is a 


and astrophysical phenomena In heavy ion 
reactions near the Fermi energy, different size isotopes, 
from Z=3 to those with Z up to Z of fission fragments, 
are copiously produced and the isotopic yields, which 
are a function of temperature and density at the time of 
the isotope formation, carry information on the 
symmetry energy of the emitting source. However the 
production mechanisms of these isotopes are not 
always completely understood. Global characteristic 


features of the experimental observables, such as 


Intermediate mass fragment production, Differential chemical potential, Double ratio temperature, 


multiplicities, mass or charge distributions and energy 
spectra or the mean energies of the fragments, have 

both 
and by dynamic 


been well reproduced by statistical 


multifragmentation | models" ^! 
transport models'**!, though the two are based on 
quite different assumptions. 

In order to further elucidate the reaction 
mechanisms and the characteristic nature of the 
emitting source, detailed experimental information, 
such as temperature and density of the emitting source, 
is indispensable. Temperature has been experimentally 
extracted mainly in two methods in this energy domain. 
One utilized the relative population of excited states of 
fragments" ?!. The other method utilized the double 


12-17] 


yield ratios of isotopes! Both results were 
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consistent near the Fermi energy. However when the 
incident energy increases, these two methods can give 
different results. In the study of Kr+Nb at 
35-1204 MeV by Xi et al, the relative population 
temperature shows more or less a constant temperature 
of around 4 MeV, independent of the incident energies, 
whereas the double ratio temperatures, derived from 
He and Li isotopes, steadily increase from 4.7 to 
9 MeV as the incident energy increases. This increase, 
however, is not observed for the double ratio 
temperature of Li and C isotopes, which shows a 
rather constant temperature of around 5 MeV. Similar 
results were reported in Aut+Au reactions at 50-2004 
MeVU*"l From these Observations, Trautmann ef 
al! concluded that two thermometers probe the 
temperature at different stages of the reaction. 

18-21] the density of the 
emitting source was also extracted, together with the 


In our recent works! 


temperature, utilizing a coalescence model technique. 
The model was applied to the light charged particles 
(LPC's) from the reactions of Part Sn and 
9 7p4l > 4Sn at 474 MeVU??8 A 4r detector, 
NIMROD, was used in the experiment?” The 
temperature of the emitting source was evaluated from 
the double yield ratios of Y(d)/Y(t) and YCHe)/Y(a). 
The density was extracted, employing the thermal 
*81 The extracted 


temperature and the density show more or less a linear 


coalescence model of Mekjian! 


correlation, that is, as the temperature decreases from 
11 MeV to 5 MeV, the density decreases from 0.29 to 
0.03po, where po is the normal nuclear density. 

The low density symmetry energy coefficients 
extracted were quite different from those extracted in 
previous works in a similar incident energy range, 
Esym~15-20 MeV and p/po-0.3-0.6"?l. Therefore it is 
very interesting to extend the previous study to heavier 
isotopes. The IMF energy is significantly modified by 
the Coulomb repulsion and possible radial flow. The 
yield of IMF's is also modified by sequential 
secondary decay process as discussed in our recent 


30-32) Therefore in order to study the density and 


works! 
symmetry energy coefficient for IMF emission, we 
derived, from the data, theratio of the difference in 
chemical potential between neutrons and protons to 
the temperature, (15745) T, which is rather insensitive 


to the secondary cooling process. The differential 


chemical potential, (j4—45), is also closely related to 
the symmetry energy at the time of the fragment 
formation, i.e., (u,-4,)-4E,.(p)ó B3] as discussed 
below, where 5=(p,-p,)/(e,+p,). Using our data 
together with results of calculations employing the 
quantum statistical model (QSM) of Hahn and 
StóckerP^, — the density, 
temperature and symmetry energy coefficient are 
studied. 


correlations between 


2 Experiment 


K-500 
superconducting cyclotron facility at Texas A&M 


The experiment was performed at the 


University. 67071 and “Ni beams were used to 
irradiate SANI He 197 Au and ??Th targets at 
404 MeV. A detailed description of the experimental 
setup and event classification has been given in 
Refs.[30,31]. Intermediate mass fragments (IMF's) 
were detected by a Si telescope placed at 20°, which 
consists of four quadrant Si detectors. Typically 6—8 
isotopes for atomic numbers, Z, up to Z=18 were 
clearly identified with the energy threshold of 
4-10A MeV, using the AE-E technique for any two 
consecutive detectors. The multiplicity of neutron and 
light charged particles (LCP) associated with each 
IMF were also measured by 16 DEMON detectors and 
16 CsI(TI) detectors around the target. 


3 Results 


3.1 Moving source fit 


The yield of each isotope was evaluated, using a 
moving source fit, where a single intermediate velocity 
(NN) source with a smeared source velocity around 
half the beam velocity was used, coupled with the 
AMD simulations as a reference. A detailed procedure 
of the multiplicity extraction was given in Ref. [31]. 
For neutrons and LCPs, three moving source fit was 
applied, i.e., projectile-like source, NN source and 
target-like source. In the following analysis the 
multiplicity of the NN source is used for all particles as 
their multiplicity. 

As direct outputs of the moving source fit, the 
source size A,=(N,+Z,) and asymmetry parameter, 
ÓNN-(N,—Z,yA, are evaluated from the extracted 
multiplicity as N,-XM;N; and Z,-XM4j, where the 
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summation is taken over all measured particles, 
including neutrons, LCPs and IMF with Zup to 18. N; 
and Z; are neutron and proton number of the particle i. 
In Fig.l, the extracted A, and óyy are plotted as a 
function of Zmax, the upper limit of the summation. A, 
and óyy values at Zmax=2, for example, indicate the 
source size and asymmetry when only light particles 
are taken into account. The emitting source size of the 
NN component varies from 45 to 80 for the small (Ni) 
to the large mass (Th) target. This indicates that all NN 
source particles originate from a rather small source, 
about a half of the projectile for the Ni targets to a 


slightly large size of the projectile for the heavy targets. 


The observation of a rather small size of the emitting 
source is consistent with previous observations, which 


indicates that, as a function of temperature, 


experimental observables such as the isoscaling 
parameter, symmetry energy coefficient and pairing 
coefficient, reflect similar reaction mechanisms , even 


though the target sizes are quite different in the 


reactions studied" 7?! 
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Fig.1 (Color online) Source size (top) and asymmetry dyy 
(bottom) of the NN source for reactions with the “Ni beam. 
Different targets are shown by different colors as indicated. 
Horizontal bars on the bottom right indicates the calculated dyy 
values from the source of 1 to 1 mixing of projectile and target 
nucleons and from the whole given reaction system, indicated 
by 1:1 and sys, respectively. 


Within errors, the evaluated asymmetry 
parameter, yy, is consistent with the values calculated 
assuming either from 1 to 1 mixing of the projectile 
and target nucleons or mixture of all nucleons. 
However the error bars, mainly originating from the 
moving source parameterization, are too large to 
determine the specific experimental asymmetry value 
for a given reaction system. Therefore in this work the 
calculated values of the 1 to 1 mixing source are used 
as Oyn for the following analysis. 


3.2 The ratio (u,-up)/ T 


The ratio of the difference in nucleon chemical 
potential relative to the temperature, (15,745)/T, can be 
experimentally determined from the isotopic yield 
ratios, using the Modified Fisher ModelP??"!. In the 
Modified Fisher model, the fragment yield of A 
nucleons with -N—-Z, Y(A, I), is given by 


Y(4,I) - CA* exp{[(W (4,1) u,N + u,Z)/T] 


l 
+ NIn(N/A) + Zln(Z/ A), M 


where C is a constant. The A” term originates from the 
entropy of the fragment and the last two terms are 
from the entropy contributions for the mixing of two 
substances in the Fisher Droplet Model!**!. W(A,D) is 
the free energy of the cluster at temperature 7. In our 
application W(A,l) is 
generalized | Weiszàcker-Beth 


formulaP?49] 


given by the following 


semi-classical mass 


When the yield ratio of two isotopes between 
=I and FH2 with same 4 is 
R(I*2,1,A) 2 Y(A,1 2) Y(A,I) is given by 

R(I 2,1, A) - exp(W (142, 4) -W (1, A) 
*(u,-u,)MT-S (142.4) D 
- Sax (1,4) 


taken, 


where S,,,(1,4)=NIn(N/A)+ZIn(Z/ A). From these 
relations, one can get 
R(I +2,1, A) - exp([u, — 4, + 2a, (Z -1)/ A^? 
-4Aa, (141) A-ó(N*1,Z-1) (3) 
-é(N,Z)/T - A(I1*2,1, A}, 


where A( *14)- S, 


the Coulomb and symmetry coefficient in W(,4) and 
Ó(N, Z) represent the pairing term. 


ü + A)S nix 


(IA), ac and asym are 
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When we focus on the yields of isobars with 
I——] and 1, one can get 


In[R(L-L 4)] - |(u, - 4) *2a.(Z -1)/ 4^ |/T. (4) 


This equation was initially applied to the € 7n "Sn 
reaction. In[R(1,-1,4)] values were calculated from the 
multiplicities of isotopes in the atomic number range 
of 3<Z<18. (u-up)T-0.71 and a,/T=0.35 were 
extracted? "!, Assuming the same temperature and 
density at the time of the fragment formation for 
different reaction systems, which results in the same 
Coulomb energy contribution in Eq.(4), (u;-up)T is 
evaluated as: 


(4, -4,)/T -(In[R,-,4)) -Qa.(Z -0/ 4^ /T).. (5) 


Here (In[R(L-1,4)] and (2a,(Z-1)/ 4'"/T) are 
the averaged experimental value over a certain range 
of Z of isotopes. Isotopes with 3€ Z <12 are used, in 
order to minimize the errors between the reaction 
Systems, since the statistics of higher Z isotopes 
becomes low for some of the reaction systems in the 
experiment. The resultant (u,—,4,)/T is shown in 
Fig.2 as a function of the asymmetric parameter 
Oy =(N-Z)/A_ of the emitting source. The values 
are well scaled by Onn. (u, — 4,)/ T 21.1ó,, +0.21 
is obtained as shown in the figure. 
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Fig.2 (Color online) (u,-u5)/T vs dyy from all 13 reactions. 
Red, black and blue represent for “Ni, Zn and “Zn beams, 
respectively. Results for different targets are also indicated by 
different symbols as shown. Data are fit by a linear function and 
the parameters are shown in the figure. 


3.3 Temperature 


In order to determine the temperature at the time of the 


double 
thermometer between different IMF isotopes was used. 


fragment formation, the isotope ratio 


The double isotope ratio temperature is generally 
given by 

T - B/ln(aR), (6) 
where R is the double yield ratio of four isotopes, B is 


a binding energy parameter and a relates to their 


[20 


spins""!, There are many choices of the double ratio 


thermometers. Here we examined the following 


thermometers"), 


13.3 MeV 
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Fig.3 (Color online) Primary temperature vs emitting source 
density. Symbols are from the double ratio temperatures and 
typical errors are given only for the results from Tana; The 
others have a similar errors in each point. Solid line is the result 
from (Hr Hp)/T fit. 


The calculated values from the different thermometers 
represent the temperatures of the final products, which 
may be significantly modified by the secondary 
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statistical decay in the cooling process. The secondary 
decay effects are carefully examined, using the QSM 
calculation as described in Ref.[25]. The resultant 
primary temperatures extracted are fitted by a 
polynomial function of the secondary temperature for 
a given density. In Fig.3 the calculated primary 
temperature for these different thermometers are 
shown as a function of the density of the emitting 
source. The secondary temperatures used are those 
averaged over all 13 reactions. For all cases the 


primary temperature increases when the source density 


increases. 
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Fig.4 (Color online) (bottom) Experimental and calculated 
(4i-Hp)/T values as a function of dyy. T=5.2 MeV is used in the 
figure. The different symbols represent the values for different 
given densities as indicated. (top) Chi-square values as a 
function of the density. 

For a given temperature and density in Fig.3, 
the differential chemical potential, (4, — 4,)/T , is 
calculated, using QSM and compared with the 
experimental values in Fig.2. An example is shown in 
Fig.4 for the case 7=5.2 MeV. For a given density and 
Onn, the differential chemical potential, (4, — u,)/T , 
is calculated. The density is varied from 0.050, to 
0.7po in this case. For a given density, (4, —4,)/T 
increases linearly as dyy increases, as shown in the 
bottom of Fig.4. From the chi-square fit, p—0.3po is 
extracted for the case. 7«4.5 MeV and 756.0 MeV, no 
solution is found to get the minimum chi-square value 


below 0.4. T and p/po values at the minimum chi 
square are plotted by a solid curve in Fig.3. From the 
(4,7 4,)/T 
temperature range of 4.5<7<6.0 MeV was evaluated, 


results of analysis, the possible 
but no specific density range was able to be 


determined in these methods. 


3.4 Symmetry energy 


The energy per nucleon of an asymmetric nuclear 
matter is approximately given by 

E(p,6)=E(p,0)+ E, (p)ó". (13) 
On the other hand the differential chemical potential is 
given by 

(1/2)(u, - u, ) - 8E(p.6)/ 06. (14) 
From these equations, one can get 
(4, — Hy) = AE (P) D3 In Fig.2, the experimental 
data are fit by a linear function and 
(4, —#,)/T =(11.141.4)6-0.2 is obtained. Using 
(4,-4,)-4E,.(p)ó from Eq.(14) and 4.5<T<6.5 
MeV, Esym coefficient is evaluated as Esym=14.6+3.5 
MeV. 


4 Conclusion 


The differential chemical potential relative to the 
temperature, (u,—4,)/T and the initial temperature 
T are studied, using the QSM calculation. The 
experimentally extracted (54,—,4,)/T values from 
the reactions studied show a linear function of Oxy 
and (u4,-4,)/T =(11.1+1.4)6,, —0.2 was extracted. 
By comparison to the results of QSM calculations, 
possible 7-p values for the emitting source were 
extracted. The extracted values are consistent with 
those of the initial temperature from the double isotope 
temperature of d, t, h and a isotopes. The extracted 
temperature range is 4.5<7<6.0 MeV. The symmetry 
energy coefficient E, is also evaluated from the 
slope value of 11.1+1.4 in the extracted experimental 
relation (uw, - 4,)/T =(11.141.4)6,, —0.2 and in the 
above temperature range. Esym=14.6+3.5 MeV was 
obtained. 
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